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Delocalized-localized transition in a semiconductor two-dimensional 
honeycomb lattice 
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We report the magneto-transport properties of a two-dimensional electron gas in a modulation-doped Al- 
GaAs/GaAs heterostructure subjected to a lateral potential with honeycomb geometry. Periodic oscillations 
of the magneto-resistance and a delocalized-localized transition are shown by applying a gate voltage. We ar- 
gue that electrons in such artificial-graphene lattices offer a promising approach for the simulation of quantum 
phases dictated by Coulomb interactions. 



Magneto-transport properties of two-dimensional elec- 
tron systems confined in semiconductor heterostructures 
and subjected to quantizing in-plane potentials were ex- 
tensively investigated^ ^. Recently renewed interest on 
these systems was triggered by theoretical proposals for 
the realization of quantum simulators in scalable solid- 
state systems^J^. These authors discussed the evaluation 
of Coulomb interaction strength and hopping properties 
of electrons confined in square arrays of coupled GaAs 
quantum dots or at the minima of standing waves created 
by surface acoustic waves. Such quantum simulators can 
represent a viable strategy to investigate quantum phases 
and many-body effects in a controlled wa}'^ ^^ and may 
open the way to the emulation of quantum phenomena 
not easily accessible in nature. To date, however, experi- 
mental tests of quantum simulators were carried out only 
in few systems including in particular cold atom4^^^^ and 
trapped ions^. 

Here we report on the fabrication and magneto- 
transport characterization of a scalable quantum simula- 
tor obtained in a modulation-doped high-mobility GaAs 
semiconductor heterostructure. The structure consists 
of an array of pillars with honeycomb symmetry that 
are formed on the surface of the semiconductor by shal- 
low etching. Owing to the dependence of band-bending 
profiles on GaAs cap-layer thickness^, these nanofabri- 
cated pillars induce a lateral potential modulation (Fq) of 
amplitude of few meV's acting on the electronic system. 
This lateral potential manifests itself by the appearance 
of a modulation of the magneto-resistivity that is peri- 
odic in the magnetic field. 

We shall show that a transition from a conducting to 
an insulating regime can be induced by changing the volt- 
age of a metallic gate deposited on top of the array. Our 
experimental findings suggest that while electrons occupy 
delocalized states at low values of the gate voltage {i.e. 
when carrier density is high and Vq is small), they lo- 
calize in the confined levels of the potential minima un- 
derneath the pillars as we sweep the gate voltage and 
increase |Vb|. This interpretation is supported by theo- 
retical evaluations of the hopping parameters for these 
quantum-confined states. 
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FIG. 1. (a) Sketch of the lateral potential with relevant geo- 
metrical parameters. J a denotes the hopping parameter and 
a = {n, m} a set of quantum numbers (see text), (b) Scanning 
electron microscopy image of the nanopatterned modulation- 
doped GaAs/AlGaAs sample with honeycomb geometry. The 
inset shows a picture of the device with four Ohmic contacts. 
(c) Calculated bound-state energy levels E^ as functions of 
|Vb|. Different curves correspond to different values of the 
angular momentum m: m — (blue circles), m — ^1 (red 
triangles), m — ±2 (green hexagons), and m — ±3 (yellow 
squares). Curves corresponding to m values which differ only 
by a sign are degenerate. For a fixed value of Vb curves with 
same symbol and color coding but higher energy correspond 
to excited states with different radial quantum numbers, (d) 
The hopping integrals J^ as functions of |Vo|. Symbol coding 
is the same as in panel (c). The horizontal dashed line labels 
the thermal energy at T = 250 mK. Vertical dashed lines 
mark the threshold values of Vb at which a new bound state 
appears. The color of the dashed lines is related to the value 
of m as in panel (c). 



The demonstration that delocalized and localized elec- 
tron states can be realized and controlled in these ar- 



tificial honeycomb lattices with the same topology of 
graphene^^ opens up exciting opportunities for the study 
of unusual electronic phases and collective modes induced 
by the lattice potential at ultra-low temperatures and for 
the simulation of massless Dirac fermion^^ ground states 
in the presence of external fields. Of particular relevance 
is the possibility to explore a regime of strong correlations 
that is not easily achievable in graphene and is charac- 
terized by inter-pillar hopping being much smaller than 
the on-site Coulomb energy. 

The sample used in this study contains a two- 
dimensional electron gas (2DEG) in a 25 nm wide, 
one-side modulation-doped Alo.1Gao.9As/GaAs quantum 
weh (QW). The 2DEG is positioned 170 nm underneath 
the surface (the doping layer is at 110 nm). It has low- 
temperature electron density Uq = 1.1 x 10^^ cm~^ and 
mobility of 2.7 x 10^ cm^/(Vs). The in-plane potential 
modulation was achieved by defining an array of Nickel 
disks (with diameter 2ro) arranged in a honeycomb- 
lattice geometry (with lattice constant a) by e-beam 
nanolithography and then by etching away the material 
outside th e disk s by inductive coupled reactive ion shal- 
low etchingl^^. Different etching depths d were realized. 
Results presented here were obtained with d = 60 nm [see 
Fig. [ija)]. Larger etching depth lead to electron local- 
ization even at zero gate voltage. The resulting lateral 
nano-patterning after removal of the Nickel, which ex- 
tends on a 100 /im x 100 fim square region, is shown 
in Fig. [lib). Experimental values of the parameters are 
a ^ 130 nm and tq ~ 30 nm. 

In order to probe the transport properties of the device, 
Ni/AuGe/Ni/Au (10/180/10/100 nm) annealed ohmic 
contacts were evaporated at the corners of the array 
[see inset of Fig. [l 'b)]. Other contacts were realized off 
the array where the 60 nm-thick GaAs layer was uni- 
formly etched away. Thanks to the low level of damage 
introduced by this process^^, the 60nm-thick pillars of 
the honeycomb array act as a weak attractive potential 
(Vb < 0) for electrons. The absence of any signature of 
conductive behavior outside the array (data not shown) 
suggests that when the etching affects the whole sample 
surface it is sufficient to fully deplete the QW. We there- 
fore estimate Vq ~ E^ ^ 4 meV where E^ is the Fermi 
energy of the original unperturbed 2DEG. A 120 nm- 
thick Al gate was deposited on top of the device to allow 
the tuning of the charge density and Vq. Electrical mea- 
surements were performed by a lock-in technique in a 
filtered ^He magneto- cryost at with 250 mK base temper- 
ature. While the overall behavior was the same after dif- 
ferent cool-downs with different illumination procedures, 
the values of the charge density at a given gate voltage 
were found to change from cool-down to cool-down. 

Figure [2|a) shows a representative low-temperature 
longitudinal magneto-resistance characteristic measured 
at gate voltage Vg = —1.62 V. The observation of dis- 
tinct quantum Hall signatures with an estimated electron 
density of 4.5 x 10^^ cm~^ demonstrates the existence 
of delocalized states populated by electrons at this den- 
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FIG. 2. (a) Longitudinal resistance Rxx as a function of the 
magnetic field at a gate voltage Vg = —1.62 V and charge den- 
sity n = 4.6 X 10^° cm~^. (b) A magnified portion of data be- 
tween 1.1 T and 1.5 T. A polynomial fit of the main Rxx curve 
was subtracted to enhance the visibility of magnetic-field- 
periodic oscillation. A scatter-plot of the resistance minima 
position is also shown, (c) Device oscillation power spectrum 
obtained by applying the Fast Fourier Transform algorithm 
on a set of dRxx/dB taken at several gate voltages (from —1.2 
to —1.62 V, corresponding to n = 4.8 and 4.5 x 10^° cm~^, re- 
spectively) . Contributions from all curves have been averaged 
to reduce the spectral noise. 



sity. Weak modulations of the resistivity [highlighted in 
Fig. ^b)] in the magnetic-field range between filling fac- 
tors u = 2 and u = 1 are present and appear only below 
a certain gate voltage. In Fig. |2|c) we report the aver- 
aged Fast Fourier Transform (FFT) relative to a set of 
dRxx/dB traces taken from Vg = —1.2 V to —1.62 V. 
Two main peaks are observed at 27 T"^ (FWHM 6 T-^) 
and 67 T~^ (FWHM 20 T~^) corresponding, respectively, 
to 2.5 and 6.5 times the geometrical area of the hexag- 
onal unit cell. Following early works based on a semi- 
classical analysi^^*^, the FFT peaks can be attributed 
to the quantization of periodic orbits, although a pre- 
cise identification requires a more extended theoretical 
investigation. 

We now address the evolution of the carrier density as 
a function of l^ in a set of two- wire quantum Hall mea- 
surements. The results presented in Fig. |3] display the 
expected decrease of the charge density above a thresh- 
old gate voltage whose value is associated with the de- 
pletion of the residual charge trapped in the Si doping 
layer. We remark that as we lower Vg the resistance be- 
comes so high that it hinders a reliable estimate of the 
electron density. This is consistent with a transition to 
a regime of electrons localized in the minima of the hon- 
eycomb potential as supported by the finite density of 
2.5 — 3 X 10^^ cm~^ found in the transition region at 
l^ ^ — 1.5 V . We remark that the two- wire configura- 
tion setup allows us to estimate the charge density even 



in this transition region. This is otherwise impossible 
because of the non-hnear behavior of the voltage-probe 
contacts in that regime. 

In order to investigate in more detail the transition 
to the localized-insulating regime we present in the in- 
set to Fig. Is] the sample resistance as a function of Vg 
at zero perpendicular magnetic field. A transition as- 
sociated to an abrupt increase of the sample resistance 
leads to a persistent insulating state. While the transi- 
tion threshold was observed to depend weakly (±200 mV) 
on the cool-down/illumination procedure, the qualita- 
tive behavior remains the same from cool-down to cool- 
down. The sample can be driven again to the con- 
ductive/delocalized regime by means of a thermal cycle 
(250 mK - 30 K - 250 mK) or by illumination. This 
latter behavior can be related to the presence of several 
trapping mechanisms largely due to the device fabrica- 
tion steps. 

The insulating behavior highlighted in Fig.[3]can be ex- 
plained in terms of electron localization in the confined 
levels of the potential minima as shown in Fig. Ilia). To 
confirm this picture, which requires negligible inter-site 
hopping between the confined levels of neighboring min- 
ima, we evaluated the hopping integrals between states in 
two pillars placed at a relative center-to-center distance 
equal to a. To this end, the external potential created 
by a single pillar was modeled by a simple rotationally- 
invariant piecewise-constant function of the 2D position 
^* Kxt(^) = "1^0 if ^ = I''"! < ^0 and zero otherwise. 
The resulting eigenfunctions are particularly simple in 
polar coordinates (r, 6) , where they can be written aa^^ 
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FIG. 3. Carrier density (after a representative cooldown) as 
a function of gate voltage I4 at T = 250 mK. The inset 
shows the behavior of the longitudinal resistance measured 
via a two-terminal setup as a function of Vg at 250 mK and 
zero magnetic field. The localized regime is achieved below 
K = -L5 V. 



^a{r) = ^n,m{r)xm{0), Xm{0) = e'^'/V^TT, with m a 
relative integer, being the eigenstates of the angular mo- 
mentum operator L^ = —ihd/dO along the z direction. 

The spectrum of bound states with energy E^ < can 
be found by solving the following transcendental equa- 
tion: 






(1) 



Here k = [2m^{\Vo\- \E^\)/h'V^', n = {2m^,\E^\/h^Y/\ 
m\) is the band mass, and Jm{x),Km{x) are standard 
Bessel functions [J^(x) and K'^{x) begin their deriva- 
tives with respect to x\. The spectrum E^^ obtained by 
solving Eq. (IT]) is reported in Fig. mc) as a function of 
Vb for various values of the angular momentum m. 

The hopping parameter Jo, , which provides the rate at 
which tunneling occurs from a pillar to a neighboring one, 
can be calculated from the following two-center integral 



^ a — 



J d^r <^l{r)V,^t{r)^a{r - a) 



(2) 



Here ^^(r) are the eigenstates of a single pillar centered 
at the origin and a is a nearest-neighbor position vec- 
tor. Numerical results for J^ are reported in Fig. Illd). 
We see that in the relevant range of potential depths, 
Vq ^ 5 meV, the hopping integral associated with the 
two lowest bound state is smaller than the thermal energy 
(horizontal dashed line). In the presence of a magnetic 
field B = Bz it is easy to show via a gauge transforma- 
tion that 



^q; — 



Jd^r e'"^ Kir)Ve^t{r)^a{r - a) 



(3) 



where we have defined (p{r) = B-{rxa)/2 and cpo = hc/e 
is the fiux quantum. The complex phase factor in Eq. (Isl) , 
which resembles the Peierls phase, and the magnetic-field 
dependence of the eigenstates ^^(r) lead to a further 
(exponential) suppression of Ja (data not shown). 

In conclusion, we have studied electron transport in 
a high-mobility GaAs quantum well in which electrons 
are subjected to a lateral artificial potential with hon- 
eycomb geometry. The possibility to control the local- 
ized/delocalized nature of the electronic states suggests 
further use of this system for quantum simulations in a 
regime dominated by electron correlations. 
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